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The exact expression for the density matrix of the kink ground state of the ferromagnetic XXZ 
chain is obtained. Utilizing this, we exactly calculate various correlation functions such as the longi- 
tudinal and transverse spin-spin correlation functions, and the ferromagnetic and antiferromagnetic 
string formation probabilities. The asymptotic behaviors of these correlation functions are also an- 
alyzed. As a consequence, we find that the spin-spin correlation functions decay exponentially for 
large distances, while the string formation probabilities exhibit Gaussian decay for large strings. We 
also evaluate the entanglement entropy, which shows interesting behaviors due to the lack of the 
translational invariance of the state. 
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I. INTRODUCTION 

The effects of quantum fluctuations of interacting 
quantum systems can be investigated by studying the 
correlation functions or quantifying the entanglement of 
the system, which are currently paid much attention. 
The exact evaluation of the correlation functions, how- 
ever, is still a challenging problem even when models arc 
completely integrable. The spin-1/2 XXZ model in one- 
dimension is one of the most fundamental model, which 
can be exactly solved by the Bethe ansatz. As con- 
cerns for the correlation functions in the antiferromag- 
netic ground state, a few exact results are known so far: 
several short distance spin-spin correlation functions (see 
[l[ and references therein) and the ferromagnetic string 
formation probability (which is the probability to find a 
ferromagnetic string with certain length) for A = 1/2 
(A: anisotropy parameter) [1, Q • 

In this paper, we intensively consider the correlation 
functions for the ferromagnetic regime of the XXZ chain. 
In this regime, it has been well-known that there are two 
translationally invariant ground states up and down, the 
state with all spins up and the state with all spins down. 
In addition to these trivial ground states, two classes of 
non-translationally invariant ground states kink and an- 
tikink were found in Q (see also @ for finite XXZ chain 
with boundary magnetic field) . Though it is not obvious 
that the kink is the ground state in the infinite lattice 
limit, the authors in [J] proved it under the assumption 
that the ground states should be "frustration free" , i.e. 
minimize not only the energy of the total Hamiltonian, 
but also the energy of the local Hamiltonian. Further- 
more, it was shown that the "frustration free" ground 
states up, down, kink and antikink are the complete set 
of the ground states [1, 0] ■ In 01 , the exact value of the 
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spectral gap was obtained and shown to be independent 
of the reference ground state. 

More recently, from the interest in the correspondence 
between the ground state of the ferromagnetic XXZ chain 
and the "quantum" Hamiltonian of the crystal melting 
model, some special correlation functions such as the 
magnetization and the longitudinal spin-spin correlation 
function of the kink ground state were exactly calculated 
1- 

Utilizing the generating function developed in in 
this paper, we derive the exact expression of the den- 
sity matrix for the kink ground state of the ferromag- 
netic XXZ chain. By using this, various correlation func- 
tions can be systematically calculated for arbitrary in- 
teraction strengths and for arbitrary distances. The fol- 
lowing correlation functions are particularly calculated 
here: the transverse and the longitudinal spin-spin cor- 
relation functions, the ferromagnetic (antiferromagnetic) 
string formation probabilities which are the probability 
finding a ferromagnetic (antiferromagnetic) string in the 
kink ground state. The en tang lement entropy of the sys- 
tem is also evaluated (see [10j for finite XXZ chain with 
boundary magnetic field). Analyzing the asymptotic be- 
haviors of these correlation functions, we find that the 
the spin-spin correlation functions exponentially decay 
for the large distances. On the other hand, both the fer- 
romagnetic and the antiferromagnetic string formation 
probabilities exhibit Gaussian decay for large strings. To 
authors' knowledge, this study is the first to investigate 
systematically the correlation functions and their asymp- 
totics in the kink ground state. 

This paper is organized as follows. In the next section, 
the kink ground state of the infinite XXZ chain in the fer- 
romagnetic regime is considered. By using the generating 
function, the exact expression for the density matrix is 
derived. From this, we concretely analyze various corre- 
lation functions in section III. The asymptotic behavior 
of the correlation functions are discussed in section IV. 
Section V is devoted to conclusion. 
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II. DENSITY MATRIX OF THE KINK 
GROUND STATE 

In this section, we derive the exact expression for the 
density matrix of the infinite ferromagnetic XXZ chain 
in the kink ground state. The Hamiltonian is defined by 



H 



E 

mGZ 



[cC^m+l + ^m^m+l + A (°'m 'm+1 ~ !)] . 



where c^, a = x,y,z are the Pauli matrices acting on 
the mth site and A is the anisotropy parameter. We shall 
consider the ferromagnetic regime A > 1. Parametrizing 
A as 



A 



<? 2 + <T 



(2) 



A > 1 corresponds to < q < 1. 

It is known that there are infinitely many zero energy 
kink ground states interpolating between spin up at — oo 
and spin down at oo 0, Q • A kink ground state is the 
superposition of kinks which have the same center. The 
center of the kink is defined as the half integer-valued 
position where the number of up spins on the right of it 
is equal to the number of down spins on the left of it. 
Denote the kink ground state whose center is at j — ^ 
(j G Z) by Any \^j) can be extracted from the 

following generating function 

!*(*)>= <S> (I T)* + z- 1 q-i^ +x) \ l) x ) 

® (g) (\i)v + ^ ih+v) \])v)- (3) 



where w = z 2 and {a\q)oo '■= XXjLoO- ~ a 3 5 )- I n the 
second equality, we have used the Jacobi triple product 
identity, 

°° 2 

(q;q)oo(-wq*;q)oc(~w- 1 qi;q) oc = ^ W J gTT. (7) 

i=-oo 

Noting (*(z)|*(z)} = ET=-oc w3 (^j\^j), we have 



(8) 



(<?; q)™ 

Next we compute IIj=i -^4- appearing in the 

numerator of ©. Note that (* fc | II"=i ^xf ' = 
unless I — k — 6 where 

6 := tt{i; fe, 4) = (+, -)} - fa, 4) = (-> +)}• (9) 

Then one can see 

n n 

n 4 £j i*w> = e i n i*i*i>. 

3=1 i 3=1 

■= o. (10) 

By induction, one can show that the following holds: 



nKi) [eH E(--) fc ETr%_, ;) 

,=i fe=o ,=1 u ¥iw W 



, (11) 



l^j) is the coefficient of z- 7 of the expansion of ^(z)), i.e, 

|*(*)> =£V|*i>. (4) 



Let us calculate the form factors, 



^in" = i4 g i^) 



(5) 



where {e 3 }, {e'j} G {+, -} and Ef* = (l±o*)/2, ^ ±T = 
cr^ = ((tJ±ktJ)/2. Xj is the position of the site where the 



operator £?S 3 acts on, and is assumed to be Xj ^ Xk for 
j ^ k. First we calculate the norm (W n |$ n ) appearing ;t1 
the denominator of (J5j) . It can be obtained by calculati 

(*(*)!*(*)>: 



m 

ing 



(*(«)!*(*)) =(-iygi;g)oo(-^ ^jgjc 
1 



(<?;<?)c 



E wJ « J2 ' 



(6) 



where Cj = g^+^j and [++] = 1, [ ] = 0, [H — ] = 

[ — h] = | (see Appendix for the proof of (fTTjO . Com- 
paring (|10[) and (jTTj) . we obtain 

71 Tl OO 



J'=l 



.7=1 



fe=0 



n /-k+n — 1 



E 



9 2 ■ (12) 



-niW(Ci-o) 

Combining |5j), (fSj> and (fT2|) . we finally arrive at 

Tt Tt OO 

^(Il4 e 0i±£ = n(0) [ ^ 1 E(- 1 ) fe 



i'=l 



fc=0 



X 2^ — J — — Q 2 8 . (13) 



In particular, when the operator ri?=i 3 preserves 
the total spin, i.e. 6 — 0, and the center of the kink is at 
— 1/2, i.e. i = 0, the density matrix elements of the kink 
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(iii) Transverse spin-spin correlation function: 



FIG. 1: Magnetization {&%). 



ground state whose center is located at —1/2 is given by 



II/-' 



>4 J )o — 



/-fe+n — 1 



E(-D fe E 



fc=0 



^niw-(o-o) 



(14) 



III. CORRELATION FUNCTIONS 

Here we analyze several crucial correlation functions: 
the magnetization, the longitudinal and transverse spin- 
spin correlation functions, and the ferromagnetic (anti- 
ferromagnetic) string formation probability which is the 
probability finding a ferromagnetic (antiferromagnetic) 
string in the kink ground state. At the end of this section, 
the entanglement entropy of the system is also consid- 
ered. These correlation functions are directly calculated 
by the density matrix (fl"4|) derived in the preceding sec- 
tion. Note that the correlation functions discussed here 
are for the kink ground state whose center is located at 
— 1/2. Other cases can also be treated by using (|13p . 

Let us list the explicit expressions of these correlation 
functions. 

(i) Magnetization 0: 



«) = l-2^(-l) fe qV^+^). (15) 
fc=0 

(ii) Longitudinal spin-spin correlation function [9j: 



(«) = l + 2£(-l) fe+1 ( 



fe=0 



(16) 



5>uV 



fc=0 



q x i + 2 — q^ 

(iv) Ferromagnetic String Formation Probability: 
P f (x,n) := (E++ •••£++ ,) 



TX2 + T 



(17) 



{x+j-\)k 



(v) Antiferromagnetic String Formation Probability: 

P a (x, n) 



(18) 




n: odd 

n: even 
(19) 

FIG. [1] shows the magnetization (cr*). One sees that 
the sign of the magnetization changes between x = — 1 
and a; = 0, which corresponds to the fact that the center 
of the kink is located at x = —1/2. One can also see that 
most of the spins are aligned up for x <C — 1 and aligned 
down for x 3> 1. As decreasing of q (or equivalcntly 
as increasing of the anisotropy parameter A), the slope 
of the magnetization curve around x = increases, and 
eventually will be infinite at the Ising limit q = (A = 
oo). 

The longitudinal and the transverse spin-spin correla- 
tion functions 



C"{x lt x) := «^>-«>K) 



(20) 



and (<t+(t~) are depicted in FIG.[2]for the case x\ = —10 
and various anisotropies. Both the correlation functions 
have a peak around x = 0. This characteristic behav- 
ior reflects the fact that the sign of the magnetization 
changes around x = and almost all spins are aligned at 
\x\ 1. These correlation functions decay exponentially 
for |x| 3> 1 (see section III in detail). 

The ferromagnetic (antiferromagnetic) string forma- 
tion probability P{(x, n) (P^x, n)) is the probability that 
the spins located in the region [x, x + n — 1] form a fer- 
romagnetic (antiferromagnetic) string. These correlation 
functions are depicted in FIG [3] and FIG QJ respectively. 
From FIG. [3}b and[4]-b, one can see Pf(x,m) > Pf(x,n) 
and P a (x,m) > P a (x,n) for m < n, as expected. As 
shown in section IV, both the string formation probabili- 
ties Pf (x, n) and P&(x, n) exhibit Gaussian decay for large 
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FIG. 2: Top: The longitudinal spin-spin correlation func- 
tions C zz (xi,x) := (cr^cr^) — (cr| 1 )(crj) for xi = —10 and 
— 10 < x < 50. Bottom: The transverse correlation functions 
(cr+cr~) for x x = -10 and -10 < x < 40. 



strings n 3> 1. As the anisotropy parameter becomes 
larger, the effect of Ising interaction becomes stronger 
than that of quantum fluctuation. In the limit q — > 
(A — > oo), the spins for x < and x > are all aligned 
up and down, respectively (cf. FIG.[T|). This is reflected 
in the slope in FIG. [3]-a becoming steeper, and the peak 
in FIG. HJ-a becoming sharper, as the anisotropy param- 
eter becomes larger. 

One can also calculate the entanglement entropy 
S(x,n), the von Neumann entropy of a subsystem 
[x, x + 1, • • • , x + n — 1] . It is defined as 

S(x, n) = -trp(x, n) log 2 p(x, n), (21) 

where p(x,n) is the reduced density matrix defined by 
tracing out the degrees of freedom of the environment 
outside the subsystem [x, x + 1, • • ■ , x + n — 1]: 

p(x,n) =tr E |*o><*o| = [P:t:::e:(x,n)] e . <e , =± , (22) 



FIG. 3: The ferromagnetic string formation probability 
(FSFP) P{(x,n) for n = 4 and various anisotropics q (a), 
and for q = 0.9 and various lengths n (b). 

where 

n 

p&;:::&(x,n) = ([[E% i j _ 1 ). (23) 

Shown in FIG. [5] is the entanglement entropy. As x — > 
±00, the entanglement entropy of the kink ground state 
is asymptotically 0, which is nothing but that of the fer- 
romagnetic ground state up and down. In FIG. [5]-a, we 
observe an intriguing phenomena that the peak of the 
entanglement entropy splits into two, as decreasing the 
parameter q (or equivalently as increasing the anisotropy 
parameter A). On the other hand, for fixed q, the same 
behavior can also be observed in FIG. [SJ-b as increasing 
the length of the subchain. 



IV. ASYMPTOTICS 

In this section, the asymptotic behaviors of the cor- 
relation functions derived in the preceding section are 
analyzed. 
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FIG. 4: The ferromagnetic string formation probability 
(ASFP) Pa.(x,n) for n = 4 and various anisotropies g (a), 
and for q = 0.9 and various lengths n (b). 



FIG. 5: The entanglement entropy S(x, n) for n — 8 and 
various anisotropies g (a), and for q = 0.7 and various lengths 
n (b). 



Let us first consider the spin-spin correlation functions. 
From and US), we find 



Z2 + 1 



X2— >+QO 



O + (2 + 4^(-l) fc+1 gV+M*i+!)) ( f 2+i. (24) 



fc=0 



Thus we obtain 

(«> - «>«) ^ A"( a;i )ff e -+ 1 for ^ » 1, 

OO 

A zz (x x ) = 4^(-l) fe (l - q k )q^+ k ^ + h), (25) 
k=0 



This shows that the longitudinal spin-spin correlation 
function decays exponentially. The asymptotics of the 
transverse spin-spin correlation function (|17[) is also eval- 



uated in the same manner: 

« CT * 2 ) ~ A + -(xi)q&*>+to for x 2 > I, 

°° i 
A+-( Xl ) = J2(-l) k q^ +{k+ ^ ){xi+ ^, (26) 

k=Q 

which shows that the transverse spin-spin correlation 
function also exhibits exponential decay. 

Now, let us analyze the asymptotics of the string for- 
mation probabilities (fl8|) . Using the identity 



i 



Pf(x,n) can be rewritten as 
P i {x,n)=q nx+n \l + B{n)), 



1, 



(27) 



00 n jk 

B{n) = £(_i)* g V + C»*-i>* ^ I 

k=i j=i '■UjtjK 1 1 



(28) 
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Since 



| g (n)|<^^4)^ 

fc=l i=i lii^v 1 9 



then 



fc=l 



(29) 



fc=l 



n 1 n 

2=1 ^ J 2=1 

=.9^(ci,...,cn)n^o) [£; ' e3] 7^r ^ 

2=1 ' j=-ao 

n j 

5S(Ci,-,Cn) = II 

2=1 



(32) 



Expressing the Laurent expansion of 3™ (Ci, Cn) as 



P f (x,n) ~<7 n:E+ ™ for n> 1. 



(30) 



This means that the ferromagnetic string formation prob- 
ability shows Gaussian decay for large strings. Note here 
that similar Gaussian behaviors are also seen in the an- 
tiferromagnetic ground state [H, ITU [l~2| . 

Finally we explicitly write down the asymptotics of the 
antifcrromagnetic string formation probability (|19j) : 



P a (x,ri) ~ q — + — 

[(q s + q 8 ) 

for n>l. 



q 2 + » ) n: odd 



> (31) 



n: even 



9S(Ci,-,Cn)=E^' A i n) ^'-^»)' ( 33 ) 

2=0 

we find (&, — ,Cn) satisfies the following recursion re- 
lation, 



Ai n+1) (C 1) ... ! Cn+l)=E(- 1 ) fe_i ^l A 2- (Cl,-,Cn). 
2=0 

(34) 



To prove (fTTjl is to show that A^ (£1, Cn) is 



V. CONCLUSION 

In this paper, the density matrix in the kink ground 
state of the ferromagnetic spin-1/2 XXZ chain has been 
exactly calculated. From this expression, the longi- 
tudinal and transverse spin-spin correlation functions, 
and the ferromagnetic and the antiferromagnetic string 
formation probability for arbitrary distances and arbi- 
trary interaction strengths have been systematically cal- 
culated. Analyzing them, we find that the spin-spin cor- 
relation functions decay exponentially for large distances, 
while the string formation probabilities show Gaussian 
decay for large strings. We have also calculated the en- 
tanglement entropy and observed the change of shape 
with the increase of the anisotropy parameter or the 
length of the subchain. 
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AW(Ci,...,C«) = (-i) J £ 



2+n— 1 



(35) 



\ HijtiiCi - CO 

Let us show this by induction. It is obvious that 
holds for n = 1. Suppose it holds for n. Then from (134p . 

A^.™ (Cii Cn+i) can be calculated as follows. 

^(n+l)^ ^ Cn+i) 

k n , v j 

= (-i)^ +1 EEcr 1/ " ' 



J=0; = 1 xCn+l/ n«=i(o-Ci) 

n »n+k _ /-n—lfk+l 
Si Si Sn+1 



1=1 rii=i (0 - c 

i=£l 

-l) k v 2 

ti-n^o-co' 



(36) 
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APPENDIX: PROOF OF JTJ) 



In the last equality, we used 



" ' - for n > 2. (37) 



E 



-{ Hi - Cfe) 



Let us show (fTTj) . Using ([3]), one obtains 

n 

(*(z)in4 £) i^)) 

2 = 1 



From (|36p , we can see (|35p holds for n + 1, which means 
(|5S) holds for any n. Thus, from (1521) . ([531) and (|3l3 |) . we 
obtain (jTTJ) . 
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